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We demonstrate bound states in the first TE and TM diffraction continua (BSC) in a linear 
periodic array of dielectric spheres in air above the light cone. We classify the BSCs according to 
the symmetry specified by the azimuthal number m , the Bloch wave vector directed along the 
array, and polarization. The most simple symmetry protected TE and TM polarized BSCs have 
m — 0 and /3 = 0 and occur in a wide range of the radius of the spheres and dielectric constant. 

More complicated BSCs with m/0 and — 0 exist only for a selected radius of spheres at a fixed 
dielectric constant. We also find robust Bloch BSCs with 0,m = 0. We present also the BSCs 
embedded into two and three diffraction continua. We show that the BSCs can be easily detected 
by the collapse of Fano resonance for scattering of electromagnetic plane waves by the array. 

PACS numbers: 42.25.Fx,41.20.Jb,42.79.Dj 


I. INTRODUCTION 

The scattering of electromagnetic (EM) waves by an ensemble of dielectric spheres has a long history of research 
beginning with Mie who presented a rigorous theory for scattering by a single dielectric sphere [1]. The overwhelming 
majority of papers since the pioneering papers by Ohtaka and his coauthors have considered periodical two- 

dimensional (2D) and three-dimensional (3D) arrays [6 9]. Surprisingly, less interest was payed to scattering by a 
linear array of dielectric spheres mostly restricted to aggregates of a finite number of spheres EMI. Guiding of 
electromagnetic waves by a linear array of dielectric spheres below the diffraction limit attracted more attention. 
There were two types of consideration: finite arrays fL3H18j and infinite arrays which were studied by means of the 
coupled-dipole approximation EMI- A consummate analysis of electromagnetic waves propagating along linear 
arrays of dielectric spheres below the light cone was provided by Linton et al [24j . 

It has been widely believed that only those modes whose eigenfrequencies lie below the light cone, are confined and 
the rest of the eigenmodes have finite life times. Recently confined electromagnetic modes were shown to exist in 
various periodical arrays of (i) long cylindrical rods j25H28j . (ii) photonic crystal slabs [2914321 ], and (iii) 2D arrays of 
spheres [33|, [34[. Similarly one may expect light trapping in the one-dimensional (ID) array of spheres with the bound 
frequencies above the light cone. Such localized solutions of the Maxwell equations are known as the bound states 
in the continuum (BSC) and were first reported by von Neumann and Wigner [35| for the stationary Schrodinger 
equation with a specially chosen oscillatory potential. Nowadays, the BSCs are known to exist in various waveguide 
structures ranging from quantum dots to acoustic periodic structures EMI where they are known as the 

embedded trapped modes, to photonic crystals [3(J l45l-[5Cij. The BSCs are of immense interest in optics thanks to 
the experimental opportunity to confine light despite that outgoing waves are allowed in the surrounding medium 

mUMM- 

II. BASIC EQUATIONS FOR ELECTROMAGNETIC WAVE SCATTERING BY A LINEAR ARRAY OF 

SPHERES 

In the present paper we consider a free-standing ID infinite array of dielectric spheres in air Fig. [lj In what follows 
we refer to all length quantities in terms of the period h of the array. We formulate the scattering theory by a periodic 
array of dielectric spheres in the form of the Lippmann-Schwinger equation similar to the approach developed for a 
periodic array of dielectric cylinders [27] 


La = Tq nc . (1) 

where the matrix L accounts for both the scattering matrix of the isolated sphere and the mutual scattering events 
between the spheres, Tq nc is given by the incident wave, and the column a consists of amplitudes off 1 of the multipole 
expansion of the scattering function. 

The exact expression of the matrix L was derived by Linton et al [24] for EM guided waves on a periodic array of 
dielectric spheres. For the reader’s convenience we present the equations and notations from the above reference. We 
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FIG. 1: A periodic infinite array of dielectric spheres illuminated by a plane wave (blue short arrow). The wave can be 
transmitted and reflected to discrete diffraction continua enumerated by integers m and n in accordance with Eqs. gb and 
d2U shown by red long arrows. 


seek the solutions of the Maxwell equations, which obey the Bloch theorem 


E(r + Rj) = e i/3Rj 'E(r), H(r + TLj) = 


— A/3 Ft j 


H(r) 


with the Bloch wave vector /3 directed along the array aligned with the z-axis (see Fig. [Tj) . Here R j = je z is the 
position of the center of the j-th sphere and e z is the unit vector along the array. Scattered EM fields are expanded 
in a series over vector spherical harmonics M™ and N™ am defined in Appendix A 


E(r) = Ei e i/3Rj Ei m K Mr(r - Rj) + b^(v - Rj)], 
H(r) = -i E, e^ E„>rNp(r - Rj) + bf Mf(r - R,-)] 


( 2 ) 


In series (|2]) the first and second terms presents TE and TM spherical vector EM fields, respectively. 

In absence of an incident wave Linton et al [24| derived the homogeneous matrix equation for the amplitudes a ™ 


and bY 1 : 


7-1 


J TE,l a l 
7-1 




y —i um ( r ,m jomm _i_ m Amm\ _ n 

Z TM,ZA - l^v\ a v ^vl ) — 

where summation over v begins with max(l,m), and the so-called Lorenz-Mie coefficients are given by 

v _ 31 (kR) [rji ( k 0 r)]' r=R -ji ( k 0 R ) [rj t ( kr)]' r=R 
TE}1 hi(koR)[rji(kr)y r=R -ji(kR)[rhi(kor)]' r=R ’ 

7 _ £Ji(kR) [rjiikor^'^n-jijkoR) [rji(kr)]' r=R 

™i l hi(koR)[rji(kr)]' r=R -eji(kR)[rhi(kor)]' r=R 5 

where k = ^fek^ and e is the dielectric constant of the spheres, 


AZ m = 4tt(—1 ) m i 


; V{V+1) 


1(1 + 1) 


E 

p=\l-v>\]l-\-v-\-p=even 


( -i) p giu P G(l , m; v, -to; p)s p , 


(3) 


(4) 


(5) 


tomm _ 

— 


2 tt (-1 y 


VK l + l) v ( v + 1) D= | 


1~\~17 — 1 

E 


p=\l-v\-\-l;l-\-is-\-p=odd 


i v 1 p \l^^'H(l,m;i/,~m;p)s p . 


The coefficients 


(l - v + p+ l)(l + v - p) (v - l +p + l)(l + v + p + 2) 
9lvp ~ + 2i/(2i/ + 1) 2(i/ + 1)(2i/+1) ’ 


Q(l,m] u,p\p) — ^- v / (2Z + l)(2i,+ l)(2p + l) ^ ^ _ m E ^ ) ( 


/ ^ p 

0 0 0 


(6) 

(7) 

(8) 
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are expressed in terms of Wigner 3-j symbols, 

1 

W(Z,ra;z/, -m\p) = ^ Q s {l,m\v,-m\p) 

s=— 1 


( 9 ) 


with 


and 


= —2m\p\Q(l,m\v,—m\p — 1), 

G±i(l,m;i/,-m;p) = =F \J(y ± m)[y =F m + l)p(p - !)£(/, ra; z/, -m ± l\p— 1), (10) 


Sp = V ^ K(koj)(e W + (-lTe-W), 

3 = 1 


(ii) 


where A/ m is normalization factor given in Appendix B. 

The next step is to account for an incident plane wave which can be expanded over vector spherical harmonics 0,13 

E ff (r) = E,=i Eb[«f ro Mr*(r) +pf m Np(r)], 

H ff (r) = -i ESi E->LMr(r) + sL^W]- (12) 

Here index a stands for plane TE and TM waves. 

Plm = ~^lrri' T lm (^) 5 Qjm = (^) 7 

pfm = -iFim-Kimia), q™ = iF lm Ti m (a), (13) 

k x = — A)q sin cq k y = fco cos <a, 


c 1 _ f_i \m;Z / 47r(2Z+l)(Z—m)!~ 

l * Y (1+m)! 5 

T '™( a ) = is r^ p r(c°sa), 

^m(a) = -£P z m (cosa). (14) 

For a particular case of normal incidence = 0 and a = —7r/2 we obtain the following from Eqs. (fl4|) 


rim = -rnP^ m (0),7T Zm = — — P/ m (0). 

da 


(15) 


The general equation for the amplitudes a™ and by 1 which describe the scattering by a linear array of spheres takes 
the following form 


7 — 1 r.rn Arnm _i_ im r?rrim/\ 

J TE,l a l 2 ^v\ a v ~r °v °vl ) 




7 — 1 im (^miomm \ um Amm\ _ 

^TM,l°l l^v\ a v °vl T 0*, ^ ; — P lm 


(16) 


Here the left hand term formulates explicitly the matrix L in Eq. m and the right hand term corresponds to the 
vector of incident wave \l/ i nc in the space of vector spherical functions notified by two integers l and m and polarization 
a. 


III. THE DIFFRACTION CONTINUA OF VECTOR CYLINDRICAL MODES 

Thanks to the axial symmetry of the array we can exploit the vector cylindrical modes for description of the 
diffraction continua which are doubly degenerate in TM and TE polarizations a. The modes can be expressed 
through the scalar function ip [1] 


iP m ,n(rA,z) = H£\ Xn r)e im * +ik ^ z 


(17) 
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Then for the TE modes we have 


and for the TM modes we have 


where 


and 


E z — 0, H z — VVn,nj 

rp _ ikp 1 d'lpm,n TT _ ik z d'lp T n,n 

r v 2 r flffy i r x“ri dr 5 




Xn r d(f) 
-ik 0 


E r = 
E(h = 


Xn 

dr 5 

E 

-4 

II 

ikz 

d^rn,n 

Xn 

dr 5 

ik z 

.,2 

1 dbm,n 


ik z 1 d'lpri 
Xn r dc , 


(18) 


X„ r 


TT _ —ikp 1 dlpm , n 

5 r Xn r d(j) 

TT _ ikp dtl> m , n 

5 ^ xl dr 1 


^2 

v £,n 


-v/2 — r,2 r 2 

An — ^0 ^2 

k z ,n = P + 2nn, n = 0,±l,±2, 


(19) 

( 20 ) 

( 21 ) 


In what follows we consider the BSCs in the diffraction continua specified by two quantum numbers m and n where 
the m is the result of the axial symmetry and n is the result of translational symmetry of the infinite linear array of 
the dielectric spheres. Note that each diffraction continuum is doubly degenerate relative to the polarization a. As 
a result of the interplay between the frequency ko and the wave number k ZjTl the continua can be open (x is real) or 
closed (x is imaginary). The axial symmetry of the system substantially simplifies the consideration of BSCs since 
the azimuthal behavior is specified by the integer m only. That reduces the dimensionality of the system from the 3D 
space of variables r, 0 , and z to the 2D space of r and z. 


IV. SYMMETRY CLASSIFICATION OF BSCS 

In the previous section we presented the theory for the scattering of plane waves by a periodic array of dielectric 
spheres based on the approach by Linton et al Q. If there is no incident wave we have La = 0 whose solutions are 
bound modes of the array. There might be two kinds of the bound modes. The first type of modes have wave number 
/3 > ko and describe guided waves along the array. These solutions found by Linton et al exist in some interval of the 
material parameters of spheres, dielectric constant e or radius L, and the Bloch wave number (3 (24]. The second type 
of bound modes with /3 < ko resides above the light cone (BSCs). It is much more difficult to establish the existence 
of the second type of bound states because a tuning of material parameters is required. However there might exist 
symmetry protected BSCs which are robust with respect to the material parameters. These BSCs have been already 
considered in the linear array of infinitely long dielectric rods |25l-[27l l29l . [3l], 134 1521] . 

The axial symmetry of the array implies that the matrices A and B split into the irreducible representations of 
the azimuthal number m which therefore classifies the BSCs. Next, the discrete translational symmetry along the 
z-axis implies that the respective wave number (3 specifies the BSC. At last, additional optional symmetries arise due 
to the inversion symmetry transformation Kf(x,y,z) = f(x,y, —z) for (3 = 0 and i r. It follows from Eq. ([TT]) that 
S2/c+i = 0? an d respectively from Eqs. ([5]) and (J6]) we obtain A= 0 if / + v is odd, and B = 0 if l + v is even. 
Moreover for arbitrary f3: B$ = 0 (see Appendix B). These relations establish the selection rules for the amplitudes 
a™ and h™ which determine the allowed BSC modes listed in Table I. 


Table I. Classification of the BSCs 
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p 

Type I of BSC 

Type II of BSC 


0 

( n m um \ 
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^0 

O 

II 

o' 

(0, bf),H z =0 

0 

0 

Kfe. 0),B Z =0 

(0, b" 2k ),H z = 0 

0 

0 

(0, b u 2k+1 ),H z = 0 

( a 2/c+l’ d), E z = 0 




















5 


The cCartesian components of the vector spherical functions transform under the inversion of z as follows 

- 9) = -(-1 ) l - m M^ y (0), M£(tt - 9) = (-1 

-d) = (-1 y- m N^ x je), N%(n - 9) = -(-1 ) l ~ m N ^)- (22) 

For f3 = 0 we have 

£j - R i) = £,• - Rj), 

Ej M^(r - Rj) = (-1 ) l ~ m E j M^ z (Kt - Rj). 

Ej - R i) = —Ej - R i)’ 

Ei ^r,( r - R i) = (- 1 )' - ” 1 Ei - R i)- ( 23 ) 

Then from these equations and Eqs. m one can obtain the following symmetric properties for the cartesian 
components of the EM fields collected in Table II. 


Table II. Symmetry properties of the eigenmodes with (3 = 0. 


Type I 

Type II 

Ex,y( — z) = ( — l) m+i E Xt y(z) 

E z (-z) = (~l) m E z (z) 
H x , y {-z) = (-1 ) m H x , y (z) 
H z {-z) = (-1 ) m+1 H z (z) 

E x ,y{ — z) = (— l) m E Xt y(z) 
E z (-z) = (-l) m+1 E z (z) 
H x ,y(~z) = (-1 ) m+1 H x , y (z) 
H z {-z) = (-1 ) m H z (z) 


Tables I and II are for the symmetry classification of the bound modes in the next sections. In particular, as it is seen 
from Table I for m = 0 and /3 = 0 the type I of BSCs is the pure TE modes while the type II is the pure TM modes. 
However when m/0 the BSCs are given by superposition of TE and TM polarized modes. Nevertheless each type I 
and type II of the BSCs presents a sort of polarization because of their orthogonality to each other. 


V. SYMMETRY PROTECTED BSCS 

In this section we present numerical solutions of Eq. m for the symmetry protected BSCs with m = 0 and (3 = 0 
embedded into the first diffraction continuum with n = 0. They constitute the majority of the BSCs in the array. 
The symmetry protected BSCs are either pure TE spherical vector modes (type I in Table I) with aS^ ^ 0 and b k = 0 
or TM spherical vector modes (type II in Table I) with a k = 0 and ^ 7 ^ 0- We show that the symmetry protected 
BSCs are symmetrically mismatched to the first open continuum. 


A. BSCs with m — 0 


Below we present numerical solutions for the TE BSCs with an accuracy of 10 . 


fc 0 = 4.24, R = 0.3,e=12, a? =( 0.7563 - 0.6542* | , 6 ? = 0, Z > 1 

V 0 

/ 0 \ 

0.08576 + 0.1161* 

0 

0.588 + 0.796* 

0 

\ -0.0002 - 0.0003* / 


k 0 = 5.0115, R = 0.4,e=12, a? = 


5 U l 


= 0 ,/ > 1 . 


(24) 


(25) 


These TE BSCs embedded into the lowest n = 0 diffraction continua of both polarizations are shown in Figs. [ 2 ] (a) 
and (b). Hereinafter we plot only the real parts of EM fields. There are also the TM BSCs: 
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k 0 = 4.7504, R = 0.3, e = 15, a? = 0, 6? = 


k 0 = 6.1522, R = 0.47, e = 15, a? = 0, = 


0 

—0.6017 + 0.7988* 

0 

0.0004-0.0006* 

0 \ 
-0.8718 + 0.3926* 

0 

0.267-0.1203* 

0 

-0.0013 + 0.0006* ) 


,l> 1 , 


J > 1 . 


(26) 


(27) 


Patterns of these TM BSCs are shown in Figs. He) and[2jd). Due to the immediate vicinity of the BSC ([27]) to the 
second diffraction continuum the BSC shows a large scale of localization around the spheres (see also the Sec. VIII 
below). 

The symmetry protected TE and TM polarized BSCs have qualitatively similar field structure with respect to 
E H but are not degenerate because of different boundary conditions for E and H at the sphere surface. The TE 
polarized BSC ([24]) and the TM polarized BSC ([26]) have the dominant contribution a® while the TE BSC ([25]) and 
the TM BSC (|27|) have the noticeable contribution of a® which is reflected in complication of the EM force lines shown 
in Figs. [2jb) and Fig. [2fd). From Table II one can see why the eigenmodes ([24]) - ([27]) are protected by symmetry 
against decay into the diffraction continua m = 0, n = 0 with TE and TM polarizations. From Eqs. m and m 
we obtain that the TE/TM continuum with k z $ = 0 (/? = 0) has the only H z /E z ^ 0 independent of z. The TE 
BSC has E z = 0 and odd H z so that this type of BSCs is symmetrically mismatched to both TE and TM continua. 
Similarly, the TM BSC has odd E z and H z = 0 and is decoupled from both TE and TM continua. 

Besides the fully symmetry protected BSCs from the third row in Table I (a!]*., 0) and (0, 6^)? we found a partially 
symmetry protected TM BSC f rom the fourth row of Table I: 


k 0 = 2.934, R = 0.4805, e = 15, a? 


/ 0.6826+ 0.0332i \ 

0 

-0.7291 - 0.0354i 

0 

-0.0008 / 


= 0,i>l, 


(28) 


however the TE BSCs with (0, & 2 fc+1 ) were not revealed in our computations. The TM BSC (1281) is symmetrically 
mismatched relative only to the mn — 0 and n = 0 continuum with TM polarization. Zero coupling of this BSC with 
the TE continuum can be achieved by tuning the radius of spheres. Patterns of EM fields and EM force lines for this 
TM BSC are shown in Fig. [3] 



FIG. 3: Pattern of the TM BSC ([28]) embedded into the TE and TM diffraction continua m — 0 and n — 0. 
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B. d zm degenerate BSCs with ft = 0 


For the m = 0 case the BSC solutions can be described by purely TE or TM modes in cylindrical coordinate as 
is shown in Figs. [2] and [3] The case m ^ 0 is fundamentally different from the former case. Nevertheless the above 
described mechanism for partially symmetry protected BSCs with m = 0 can be exploited for even the case m/0. 
Obviously, the system has the time reversal symmetry which implies that such BSCs are degenerated over ± m. Let us 
start with the type I BSC with m = 1 which has the odd E z and the even H z according to Tables I and II. This BSC 
is symmetrically mismatched with the TM diffraction continuum m = 1 and n = 0 which is independent of z. The 
coupling with the TE continuum can be canceled by tuning the radius. The result of computation of this partially 
symmetry protected type I BSC (a2fe?&2fe+i) is the following 


m = 1, k 0 = 2.847, R = 0.3945, {aj, b}) 


0 0.6662 + 0.4273* 

—0.33 + 0.5145* 0 

0 -0.0048 - 0.0031* 

\0 0 


l > 1 


(29) 


and is shown in Fig. H] (a). The type II BSC (« 2 fc+i• & 2 fc) with m = 2 has even E z and odd H z . It is symmetry 
protected against decay into the TE continuum with m = 2 and n = 0 and coupling with the TM continuum is 
canceled by tuning the radius with the following result 


m = 2, fc 0 = 3.086, R = 0.471, (a 2 l ,b 2 l ) = 


t 0 

—0.2142 + 0.69642 
0 
0 

V 0 


0.6545 + 0.2013i 
0 

-0.0057- 0.00182 
0 
0 


, / > 2 . 


(30) 


All components of electric and magnetic fields are nonzero and localized around the array as shown in Fig. 0| We 
show the EM field around only one sphere because the pattern is periodically repeated along the z-axis. One can see 
that the value of the azimuthal number m is reflected in the structure of force lines in the xy-plane while the number 
of the amplitudes a™ reflects in the structure of lines along the z-axis. Figure 0| clearly shows that the BSCs with 
m/0 are neither TE polarized nor TM polarized. 


VI. ROBUST BLOCH BSCS WITH /? ^ 0 AND m = 0 


Could the Bloch BSC occur at f3 ^ 0 in the continuum of free-space modes? This question was first answered 
positively by Porter and Evans [Hj who considered acoustic trapping in an array of rods of rectangular cross-section. 
Marinica et al [46| demonstrated the existence of the Bloch BSC with /3 ^ 0 in two parallel dielectric gratings and 
Ndangali and Shabanov [26] in two parallel arrays of dielectric rods. In a single array of rods positioned on the surface 
of bulk 2d photonic crystal multiple BSCs with /3 > 0 were considered by Hsu et al (29!]. The Bloch BSCs in a single 
array of cylindrical dielectric rods in air were also reported in Ref. na. Such traveling wave Bloch BSCs with the 
eigenfrequencies above the light cone are interesting because the array serves as a waveguide although only for fixed 
P (see summary of BSCs in Fig. [8]) in contrast to the bound states below the light cone [24j . 

According to Table I the Bloch BSCs with (5/0m = 0 have only the nonzero components a® or b®. Let us first 
consider type I BSCs with b® = 0 which have E z — 0 and, therefore, are decoupled with the TM continuum but 
coupled with the TE n = 0 and m — 0 continuum. We show numerically that this coupling can be canceled under 
variation of /3. The numerical results are collected in Eq. (1311) with the pattern of EM fields shown in Fig. [5j 


k 0 


3.6505, 


R = 0.4,e = 15, P = 1.2074, 


( 4 , bi) = 


( 0.1053 - 0.06382 0 \ 
0.1918 + 0.31612 0 
0.6046 + 0.55722 0 
0.7873 + 0.47772 0 
-0.0033 - 0.0054* 0 j 


l > 1 . 


(31) 


Although this BSC occurs at the fixed value of f3 there is no necessity to tune the material parameters of the spheres 
and therefore the BSC can be referred to as robust which is attractive from an experimental viewpoint. We managed 
to find only type I BSCs for e = 15 and none of type II. Such a difference between the types is related to different 
boundary conditions for electric and magnetic fields at material interfaces. 
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FIG. 4: (Color online). BSC with /3 = 0 embedded into the TE and TM continua n — 0 and m/0: (a) the Ttype I BSC (l29|) 
with m — ±1 and (b) the type II BSC (l30l) with m = d=2. 




VII. THE BOUND STATES EMBEDDED INTO TWO AND THREE DIFFRACTION CONTINUA 

According to Sec. Ill the continua in the form of outgoing cylindrical waves are specified by two numbers m 
and n which define k Z:Tl . Above we presented numerous BSCs embedded into the first diffraction continuum with 
n = 0. However there might be BSCs embedded into a few continua as it was shown for the case of grating structures 
BUS- Let us consider a TE BSC with m = 0 and /3 = tt with the nonzero components a^ k ^ 0, E z =0 and odd 
component H z according to Table I. This BSC is coupled with the TE polarized radiation continua m = 0,n = 0 and 
m = 0,n = — 1 which have k Z: o = n and k z -\ = —it respectively. Because of degeneracy of the continua we can form 
linear combinations with both even and odd H z . Then, obviously, the the BSC remains coupled with the odd H z 
continuum. This coupling can be canceled by variation of the sphere radius to give the following result for the BSC 
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FIG. 5: (M field configurations of the TE Bloch BSC with ft = 1.2074 given by Eq. (13T1) embedded into TE and TM continua 

m = 0, n = 0. 


amplitudes 


TO = 0, (3 = tt. k 0 = 5.0185, R = 0.35456, e = 15, a° 




0.1552-0.0123* 

0 

-0.9847 + 0.0776* 
0 


l > 1 . 


(32) 


The EM field and force lines are shown in Fig. El 





FIG. 6: Pattern of the TE BSC with m — 0 and [3 = n embedded into four TE and TM continua with m — 0, n = 0 and 
m = 0, n = —1 respectively given by Eq. (1321) . 

We also found a type II BSC (0, b 2 k) with m = 0 and /3 = 0 with odd E z and H z = 0 embedded into three continua 
with n = 0 and n = ±1 shown in Fig. [71 As is shown above this BSC is completely decoupled from the TE and TM 
radiation continua with n = 0 due to the symmetry. As for the other continua with n = ± 1 the BSC is decoupled from 
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FIG. 7: Pattern of the TM BSC with m — 0, f3 = 0 embedded into three continua with m — 0, n = 0 and m — 0 and n = ±1 
given by Eqs. (l33l) . 


the TE continua. Similar to the previous case the degenerate TM continua have k Zi ±i = E2tt and can be superposed 
into the continua with either even or odd E z . Thus, the Type II BSC is decoupled with the continuum with even 
E z . By variation of the sphere radius we achieved zero coupling with the continuum with odd E z with the following 
solution: 


m = 0, /3 = 0, fc 0 = 8.9129, R = 0.4274, e = 15, 6? 


( 0 \ 

-0.2273 - 0.1508i 

0 

—0.8015 — 0.532i 

0 

V 0.0082 + 0.0055i 


l > 1 . 


(33) 


We collected the BSC frequencies fco and Bloch vectors f3 in Fig. [51 


VIII. EMERGENCE OF THE BSC IN SCATTERING 

Scattering of plane waves by periodic 2D arrays of dielectric spheres was first considered in the pioneering papers 
by Ohtaka et al M (see also Ref. 6]). Scattering by aggregates of spheres was considered in the framework of 
multi-sphere Mie scattering 0Eo,I2, nevertheless to our knowledge the scattering by a ID infinite array of dielectric 
spheres has not been considered so far. In this section we present the results of numerical computations for differential 
and total cross-sections of the infinite array with the focus on resonant traces of the BSCs similar to the scattering by 
an array of dielectric rods [13, HH, <3lj]. In what follows we restrict ourselves to the BSCs which are standing localized 
solutions with /? = 0. The general theory of scattering in terms (a] 71 , bf 1 ) is formulated in the form of Eq. m which 
allows to find the amplitudes. After the amplitudes are found from Eq. m one can expand EM fields ([2]) over vector 
cylindrical modes to calculate the cross-sections. 

While the BSCs are given by the homogeneous part of Eq. (pQ) with T^ nc = 0, the scattering fields are given by the 
solution of inhomogeneous Eq. © with an incident plane wave at the right hand part. As it follows from Eqs. CT 
and m only one diffraction channel n = 0 is open for low frequencies fco where the majority of the BSCs occur. 
Although the BSCs can not be probed directly by an incident wave they are seen as collapses of the Fano resonance 
when the BSC point is approached in the parametric space. That phenomenon was observed for the scattering of EM 
waves by arrays of rods [2514271. i29l , EkO , l33i . Iddl . [49j and layered sphere [55[ . In this section we report a similar Fano 
resonance collapse in the differential and total cross-sections vs frequency when the wave number k z tends to zero or 
the radius of the spheres approaches the BSC radius. The Fano resonance for the present system can be interpreted 
as an interference of the optical paths through and between the spheres. We restrict ourselves to the BSC effects on 
the cross-section for the fully symmetry protected BSCs and the BSCs degenerate over m = ±2. 
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( 3 / 7T 

FIG. 8: (Color online). BSC frequencies and Bloch vector f3 relative to the light line ko = (3. Dash and dash-dot lines show 
thresholds where the next continua n = ±1 and n — — 2 are opened. Fully symmetry protected BSCs El and El are marked 
by open circles (TE BSCs are in red color) and TM BSCs are in green color), TM BSC (l28l) is marked by star, two degenerate 
BSCs (l29l) and (l30l) are marked by points, Bloch BSC (I3T1 ) with /3 ^ 0 is marked by a rhombus, and BSCs (l32l) and (l33l) 
embedded into two and three continua are marked by crosses. 


Let us consider an incident plane wave with the wave vector in the x, z plane and polarizations: (a) TE polarized 
with the electric field along the y-axis and (b) TM polarized with the magnetic field along the y-axis. For m = 0 and 
k z 7 ^ 0 Eqs. (fl3l) and (fl4|) give that pf 0 E = 0 and qf 0 M = 0. Then taking into account that = 0 (see Appendix B) 
we have the following from Eqs. m for the TE incident plane wave 

= (^) 

We do not present here the sectors of wave scattering with m/0 since the type I BSC belongs to the sector mm 0. 
The second equation gives b® = 0, so that the plane wave with TE polarization after the scattering is given by a® only. 
Then the type I BSCs are quasi BSCs weakly coupled with the TE continuum for small k z . That results in a sharp 
resonant contribution in the cross-section <jte,te as shown in Fig. [9] (a). The cross-sections <Jte,tm,&tm,tm and 
&tm,te have no features related to these BSCs and are not shown in Fig. [9] (a). If the plane wave is incident onto 
the array normally a = —7r/2 (k z = 0) we have a fully invisible type I BSC that is shown by dash line in Fig. [9] (a). 
Alternatively, the symmetry protected type II BSCs with the only nonzero bk can be observed via the cross-section 
&tm,tm as shown in Fig. [9] (b). Thus, although the BSCs have no effect for the normal incidence they are detected 
by the collapse of Fano resonances in total cross-sections for k z 0. 

Next, consider the effect of the BSCs with m = 2 given by Eq. (| 30 |) on the cross-section. We begin with the TE 
plane waves incident onto the array normally {k z = 0). Then we have from Eqs. (fl3|) - ([l5|) that pf 2 E = 0, q[ 2 E ^ 0 for 
odd Z, and pJ 2 E ^ 0, q[ 2 E = 0 for even Z. Therefore as Eq. (fl6|) shows there are only type II solutions for scattered waves 
with the amplitudes (a 2 /c+i, &2fe)- Table I shows that they belong to the same type of BSCs with m = 2. Therefore 
in the vicinity of Rbsc = 0.471 this BSC is coupled with the TE continuum and gives the resonant contribution in 
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FIG. 9: Total cross-section for scattering of the plane wave incident by the angle 0 onto the array, (a) Scattering of the TE 
plane wave is strongly affected by the presence of the symmetry protected type I BSC (HU) with the eigenfrequency ko = 4.24 
for R = 0.3 and e = 12. (b) Scattering of the TM plane wave is strongly affected by the presence of the symmetry protected 
type II BSC (l26l) with the eigenfrequency ko = 4.7504 for R = 0.3 and e = 15. 


the cross-section gte,te that is demonstrated in Figs. [TO] (a) and (b). As for the scattering of the TM plane waves 
there are no resonant features as shown in Fig. [TO] (b) by dash line. One can see in Fig. [TO] (c) bright features of the 
differential cross-sections near the eigenfrequency of the BSC caused by the resonant contribution of the amplitude 
A 2 at the azimuthal angles 0 = 0, ±90°, 180°: 

^7 = cro\^2A m cos(m<f>)\ 2 . (35) 

It is clear that for the sphere radius close to Rbsc = 0.471 the BSC solution dominates in the near field zone. The 
solution can be presented as 


^ = a^BSC + ^ (36) 

where a has a resonant behavior over frequency ko with the resonant width j ~ \R — Rbsc\- Analytical expression 
for the resonant width can be derived following Refs. [H, [56). Thus we have slowly decaying quasi BSC modes above 
the light cone similar to those considered in Ref. [57]. That effect is important for concentration of light by touching 
spheres [33|, IH| notified as the harvesting capability of the system. Fig. [10] (d) illustrates the harvesting capability of 
the array of spheres in the vicinity of the BSC ([30]) . Solid blue line shows the contribution of the background 0 = \ \<$>\\ 
where 11 • ■ • 11 is the norm of vector 3>. 


IX. SUMMARY 

Recently the BSCs above the light cone were shown to exist in various systems of ID arrays of dielectric rods and 
holes in a dielectric slab [25l Kill. l46j . Similar acoustic BSCs called embedded trapped Rayleigh-Bloch surface waves 
were obtained in a system of material rods jlol-Eij. One could ask why BSCs occur in periodic dielectric structures 
(gratings) but not in homogeneous structures like a slab or a rod which can support guided EM modes below the 
light cone only. Let us begin with the simplest textbook system of a dielectric slab infinitely long in the x, y plane 
with the dielectric constant e > 1. The Maxwell equations can be solved by separation of variables for scalar function 
0(x, y, z) — e lkxX+ i k yyjjAz) to result in bound states below the light cone fcg = k^ + ky [59| while all solutions above the 
light cone are leaky [60|. The situation can be cardinally changed by replacing the continual translational symmetry 
by the discrete symmetry e(x,y,z) = e(x + ph,y,z) where p = 0,=bl,±2 ,..and h is the period of the structure. 
Then the radiation continua of plane waves e lk x,nx+ik y y+ik z z are q Uan ti ze d k Xjn = /3 + 2e n/h, n = 0, d=l,±2,... 
with the frequency fcg = k l + ky + k‘j. Here /3 is the Bloch wave vector along the x-axis, and the integer n refers 
to the diffraction continua [26|. The physical interpretation of this statement is related to the slab with the discrete 
translational symmetry being considered as a ID diffraction lattice in the x-direction. Let us take for simplicity 0 = 0 
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FIG. 10: The effect of the BSC (l30|) with m = 2, ko = 3.086 and R = 0.471 in: (a) differential cross-section vs frequency 
and the azimuthal angle, (b) total cross-sections for different radii of the spheres close to the BSC radius (1291) for plane wave 
illuminating the array normally, (c) Frequency behavior of the amplitudes A m in the expansion (1351) . (d) Harvesting capability 
of the quasi BSC at R = 0.473. The dashed red line shows the contribution of the BSC into the scattering function; the blue 
solid line shows the background 0. 


and k y = 0. Assume there is a bound solution with the eigenfrequency k^^sc > 0 which is coupled with all diffraction 
continua enumerated by n. Let ko,BSC < 27r/ft, i.e., the BSC resides in the first diffraction continua but below the 
others. Because of the symmetry or by variation of the material parameters of the modulated slab we can achieve that 
the coupling of the solution with the first diffraction continuum equals zero [ 27 HH. However the solution is coupled 
with evanescent continua n = 1,2,... giving rise to exponential decay of the bound solution over the z- axis. The 
length of localization is given by L ~ 1 2 Therefore the evanescent diffraction continua play a principal 

V 4^ /h ~k Q bsc 

role in the space configuration of the BSCs. Moreover, one can see from Fig. [8] that in the limit h oc the BSCs 
with frequency kossc —^ 0 leave no room for the BSCs with k^^sc > 0. 

In the present paper we choose another strategy to quantize the radiation continuum. We replace the rod with 
continual translational symmetry by a periodic array of dielectric spheres. Because of the axial symmetry of the 
array aligned along the z-axis the quantized continua are specified by two integers, m and n. The first integer is the 
azimuthal quantum number and the second number defines discrete directions of outgoing cylindrical waves (ED given 
by the wave vector k z , n = /3 + 2iTn/h in each sector m where /3 is the Bloch vector along the array. The bottoms of 
the particular continua with m = 0 and n = 0, ±1 and n = — 2 are shown in Fig. [51 By arguments similar to those 
presented above for the grated slab we obtain that the BSC with /3 = 0 embedded into the first radiation continuum 
m = 0, n = 0 is localized around the array with the radius of localization given by 


^ Att 2 / h 2 — k\ 


2 - 

0 ,BSC 
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The symmetry of the system is also important for classifications of the BSCs which are labeled by the azimuthal 
number of the continuum m of cylindrical vectorial waves and the Bloch wave vector /3. 

(1) The symmetry protected BSCs constitute the vast majority of BSCs which are symmetrically mismatched with 
the first diffraction continuum with m = 0 and n = 0 of both polarizations. The EM field configurations of such BSCs 
presented in Fig. [2] show hybridizations of a few orbital numbers l = 2,4, 6,... which specify the BSCs as multipoles 
of high order. Therefore the BSC solutions can not be obtained by the use of the dipole approximation PH. The 
most remarkable property from an experimental viewpoint is the robustness of the BSCs relative to the choice of 
the material parameters of the dielectric spheres. We present in Fig. [3] an example of the BSC which is symmetry 
protected relative to the TM diffraction continuum but has a zero coupling to the TE continuum obtained through 
variation of the sphere radius. 

(2) By tuning of the radius of the spheres we found BSCs in the next sectors of continua with m ^ 0. These 
BSCs shown in Fig. 0]are remarkable by degeneracy over the sign of the azimuthal number. Each BSC with =b m has 
opposite the Poynting vector. 

(3) We demonstrated that the BSC can be accessed not only by variation of the material parameters but also by 
variation of Bloch wave vector /3 along the array axis. Patterns of the Bloch BSCs are presented in Fig. OH 

(4) We found the trapped EM modes embedded into two diffraction continua with n — 0 and n — 1 (Fig. [6]) and 
even three continua with n = 0 and n = ± 1 (Fig. [7j). 

The symmetry properties of the BSC play a very important role since it is difficult to provide a zero coupling even 
with the lowest continua with n = 0 because of the degeneracy in polarization. Nevertheless the symmetry allows one 
to decouple the BSC at least with some particular continua. 

The advantage of dielectric structures is a high quality factor and a wide range of BSC wavelengths from microns 
(photonics) to centimeters (microwave) as dependent on the choice of the radius of the spheres. Although the BSCs 
exist only in selected points in the parametric space there is a nearest vicinity of the BSC point where the BSC 
predominantly contributes into the cross-section and the EM field in the near field zone as seen from Figs. l9landfT0l 
That leads to extremely efficient light harvesting capabilities [58j. The far zone EM fields can also show abundant 
features related to the BSCs. In particular Fig. [8] (a) demonstrates the effect of antenna when the BSC with azimuthal 
number m = 2 converts the EM energy into the perpendicular directions. 
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X. APPENDIX A 


The solution of the Maxwell equations inside and outside of the dielectric sphere can be written via the scalar 
function = V ; ( r )^ / Zm(^,0) where the radial solution is 


,,x f cji(Vekor) \ir<R 
\ aji(k 0 r) +bh\ 1 \k 0 r) if r > R 


ji and /i[ 1J are spherical Bessel and Hankel functions defined as 


jl{x) = \ifx Jl+1/2{x) ’ h ^ {x) = \[y x JH ^ {x) - 


Yim are the spherical functions given by 

Y l m (e,<p) = (-l) m A/ m P™(cos6>)e im ^ 
Pr(x) = (1 - _ i)» 

Yr’(6^) = (-l) m Y- m (6,4>), 

and 


(37) 

(38) 

(39) 


A/m — 


I (21 + l)(l — to )! 


(40) 


47r(Z + to)! 

Following Stratton [1] we introduce two independent vectorial fields expressed through a single scalar function ip which 
satisfies the wave equation as follows 


MJ” = V x (ripi m ), NJ” = -VxMr. 


(41) 
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Then for TE vector spherical modes we have 


and for the TM vector spherical modes 



M r A 

NT* J 



(42) 


(43) 


XI. APPENDIX B 


The value B ™ is expressed via 


U(l,0,p,0,p) = Q + + Q- 

for l + v +p odd according to Eqs. (l6l)- (fTT)l) where 


Q± = T y/u(v + l)p{jp - 1)£/(Z, 0, u,H,p — 1) 

5(1,o,±i,p-1)^-#TiK^TI)M( J (J v 0 % 1 

according to Eq. it. Using the property of 3j-symbols 

f dl ^*2 JS A _ /_^\Ji+J l 11 2+J l 3 d2 JS \ 

m 2 ' V _mi -m2 -m 3 / 

we obtain 

f l up — 1 \ _ / Z ^ p — 1 \ 

U 1 - 1 J v° - 1 1 J 


(44) 


(45) 


(46) 


(47) 


if / + v +p — 1 is even. Therefore we have from Eqs. (ITU) and (|T5j) that 1~L(l , 0, z/, 0,p) = 0 and respectively, = 0. 


[1] J. A. Stratton Electromagnetic Theory (McGraw-Hill, New York, 1941). 

[2] K. Ohtaka, Phys. Rev. B19, 5057 (1979). 

[3] K. Ohtaka, J. Phys. C: Solis State Phys.13, 667 (1980). 

[4] M. Inoue, K. Ohtaka, and S. Yanagawa, Phys. Rev. B25, 689 (1982). 

[5] H. Miyazaki and K. Ohtaka, Phys. Rev. B58, 6920 (1998). 

[6] A. Modinos, Physica 141 A, 575 (1987). 

[7] J.H. Bruning and Y.T. Lo, IEEE Trans. Antennas and Propagation, AP-19, 378 (1971). 

[8] F.J. Garcia de Abajo, Rev. Mod. Phys. 79, 1267 (2007). 

[9] K.X. Wang, Z. Yu, S. Sandhu, V. Liu, and S. Fan, Optica 1, 388 (2015). 

[10] K.A. Fuller and G.W. Kattawar, Opt. Lett. 13, 90 (1988) ibid 13, 1063 (1988). 

[11] A.-K. Hamid, I.R. Ciric, and M. Hamid, Canad. J. Phys. 68, 1157 (1990). 

[12] Yu-lin Xu, Appl. Optics 34, 4573-4588 (1995) 

[13] D.W. Mackowski, J. Optics, 11, 2851 (1994). 

[14] M. Quinten, A. Leitner, J.R. Krenn, and F.R. Aussenegg, Opt. Lett. 23, 1331 (1998). 

[15] A. Quirantes, F. Arroyo, and J. Quirantes-Ros, J. Colloid and Interface Sc. 240, 78 (2001). 

[16] Pi-Gang Luan, Kao-Der Chang, Opt. Express, 14, 3263 (2006). 

[17] R. Zhao, T. Zhai, Zh. Wang, and D. Liu, J. Lightwave Tech. 27, 4544 (2009). 

[18] J. Du, S. Liu, Z. Lin, J. Zi, and S.T. Chui, Phys. Rev. A79, 051801R (2009), ibid 83, 035803 (2011). 

[19] M. Gozman, I. Polishchuk, and A. Burin, Phys. Lett. A372, 5250 (2008). 

[20] G.S. Blaustein, M.I. Gozman, O. Samoylova, I.Ya. Polishchuk, and A.L. Burin, Opt. Express, 15, 17380 (2007). 

[21] B.T. Draine and P.J. Flatau, J. Opt. Soc. Am A25, 2693 (2008). 

[22] R.S. Savelev, A.P. Slobozhanyuk, A.E. Miroshnichenko, Yu.S. Kivshar, and P.A. Belov, Phys. Rev. B89, 035435 (2014). 

[23] A. Krasnok, S. Makarov, M. Petrov, R. Savelev, P. Belov, and Yu. Kivshar, arXiv: 1503.08857vl [physics.optics] (2015). 

[24] C.M. Linton, V. Zalipaev, and I. Thompson, Wave Motion, 50, 29 (2013). 





17 


[25] S.P. Shipman and S. Venakides, Phys. Rev. E71, 026611 (2005). 

[26] R.F. Ndangali and S.V. Shabanov, J. Math. Phys. 51, 102901 (2010). 

[27] E.N. Bulgakov and A.F. Sadreev, Phys. Rev. A90, 053801 (2014). 

[28] Zhen Hu and Ya Yan Lu, J. Optics 17, 065601 (2015). 

[29] Chia Wei Hsu, et al , ”Bloch surface eigen states with in the radiation continuum”, Light: Science and Applications 2, 1 
(2013). 

[30] Yi Yang, Chao Peng, Yong Liang, Zhengbin Li, and S. Noda, Phys. Rev. Lett. 113, 037401 (2014). 

[31] Bo Zhen, Chia Wei Hsu, Ling Lu, A.D. Stone, and M. Soljacic, Phys. Rev. Lett. 113, 257401 (2014). 

[32] Chia Wei Hsu, Bo Zhen, J. Lee , Song-Liang Chua, S.G. Johnson, J.D. Joannopoulos and M. Soljacic, Nature, 499, 188 
(2013). 

[33] Mingda Zhang and Xiangdong Zhang, Scientific Rep. 5:8266, 1 (2015). 

[34] Maowen Song, Honglin Yu, Changtao Wang, Na Yao, Mingbo Pu, Jun Luo, Zuojun Zhang, and Xiangang Luo, Opt. 
Express, 23, 2895-2903 (2015). 

[35] J. von Neumann and E. Wigner, Phys. Z. 30, 465 (1929). 

[36] J.U. Nockel, Phys. Rev. B46, 15348 (1992). 

[37] C.S. Kim, A.M. Satanin, Y.S. Joe, and R.M. Cosby, Phys. Rev. B60, 10962 (1999). 

[38] O. Olendski and L. Mikhailovska, Phys. Rev.B 66, 035331 (2002). 

[39] A.F. Sadreev, E.N. Bulgakov, and I. Rotter, Phys. Rev. B 73, 235342 (2006). 

[40] R. Porter and D.V.Evans, J. Fluid Mech. 386, 233 (1999). 

[41] D.V. Porter and R. Porter, Q. J. Mech. Appl. Math. 55, 481 (2002). 

[42] R. Porter and D.V.Evans, Wave Motion 43 29 (2005). 

[43] C.M. Linton and P. Mclver, Wave Motion 45, 16 (2007). 

[44] D.J. Colquitt, R.V. Craster, T. Antonakakis, and S. Guenneau, Proc. R. Soc. A471, 1 (2015). 

http: / / rspa.royalsocietypublishing.org 

[45] E.N. Bulgakov and A.F. Sadreev, Phys. Rev. B78, 075105 (2008). 

[46] D. C. Marinica, A. G. Borisov, and S.V. Shabanov, Phys. Rev. Lett. 100, 183902 (2008). 

[47] S. Longhi, Phys. Rev. A78, 013815 (2008). 

[48] S. Longhi, J. Mod. Optics 56, 729 (2009). 

[49] E.N. Bulgakov and A.F. Sadreev, Opt. Lett. 39, 5212 (2014). 

[50] N. Rivera, Chia Wei Hsu, Bo Zhen, H. Buljan, J.D. Joannopoulos, and M. Soljacic, arXiv: 1507.0092vl (2015). 

[51] Y. Plotnik, O. Peleg, F. Dreisow, M. Heinrich, S. Nolte, A. Szameit, and M. Segev, Phys. Rev. Lett. 107, 183901 (2011). 

[52] M. Lopez-Garcfa, J.F. Galisteo-Lopez, C. Lopez, and A. Garcfa-Martm, Phys. Rev. B85, 235145 (2012). 

[53] G. Corrielli, G. Della Valle, A. Crespi, R. Osellame, and S. Longhi, Phys. Rev. Lett. Ill, 220403 (2013). 

[54] S. Weimann, Yi Xu, R. Keil, A.E. Miroshnichenko, A. Tunnermann, S. Nolte, A.A. Sukhorukov, A. Szameit, and Yu.S. 
Kivshar, Phys. Rev. Lett. Ill, 240403 (2013). 

[55] F.Monticone and A.Alu, Phys. Rev. Lett. 112, 213903 (2014). 

[56] E. N. Bulgakov , K. N. Pichugin , A. F. Sadreev , and I. Rotter, JETP Lett. 84, 430 (2006). 

[57] T. Ochiai and K. Sakoda, Phys. Rev. B63, 125107 (2001). 

[58] A. I. Fernandez-Dominguez, S.A. Maier, and J. B. Pendry, Phys. Rev. Lett. 105, 266807 (2010). 

[59] J.D. Jackson, Classiacl Electrodynamics , (N.Y. 1962). 

[60] J. Hu and C.R. Menyuk, Adv. Opt. Photonics 1, 58 (2009). 


